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Chiral Surface States in the Bulk Quantum Hall Effect 



in 

OS 



o 
O 



> 

On 
OO 
O 
O 

ON 



o 



Leon Balents and Matthew P. A. Fisher 
Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106-4030 

(February 1, 2008) 

In layered samples which exhibit a bulk quantum Hall effect (QHE) , a two-dimensional (2d) surface 
"sheath" of gapless excitations is expected. These excitations comprise a novel 2d chiral quantum 
liquid which should dominate the low temperature transport along the field (z-axis) . For the integer 
QHE, we show that localization effects are completely absent in the "sheath", giving a metallic z-axis 
conductivity. For fractional filling v — 1/3, the "sheath" is a 2d non- Fermi liquid, with incoherent 
z-axis transport and a^^ ~ T'^- Experimental implications for the Bechgaard salts are discussed. 



PACS numbers; 73.40.Hm, 75.30. Fv, 71.27.+a 

Disorder has a profound effect on transport in two- 
dimensional (2d) electron systems. |l[]. In the absence 
of an applied magnetic field, all the electronic states are 
believed to be localized due to strong quantum interfer- 
ence effects. With weak disorder, when the conductivity 
a is larger than /h, there is a weak-localization regime 
described by diffusive behavior with logarithmic temper- 
ature corrections. In the low temperature limit, however, 
a crossover to strongly localized behavior with a <C /h 
is always expected. In this regime the conductivity drops 
rapidly with temperature and vanishes as T — > 0. 

In the presence of a magnetic field, 2d localization can 
be circumvented by tuning to the center of Landau lev- 
els. At these isolated transitions between quantum Hall 
plateaus, extended states and a temperature independent 
conductivity of order /h are predicted, consistent with 
experiment [^j . There are also other 2d systems which ex- 
hibit such "metallic" behavior at isolated transitions, for 
example 2d films at the superconductor-insulator transi- 
tion. But away from such transitions, when the conduc- 
tivity is well below e^//i, it is invariably strongly temper- 
ature dependent, and insulating at zero temperature. 

In this paper, we describe a novel class of anisotropic 
2d electronic phases, which surprisingly can exhibit 
"metallic" conductivities much smaller than jh. These 
2d phases arise at the surface of bulk f/iree-dimensional 
quantum Hall samples. We consider layered samples 
which exhibit independent quantum Hall states in each 
layer when a large perpendicular magnetic field is ap- 
plied. This requires an inter-layer tunneling amplitude 
t small compared to the 2d quantum Hall gap Eg. In 
such materials, the surface of the sample is enveloped by 
a sheath of current-carrying states as depicted in Fig- 
ure 1. This chiral surface phase is the 2d analog of the 
Id states at the edges of a single layer quantum Hall 
fluid. There are currently two candidate experimental re- 
alizations of such systems: 2d electron gas multi-layers, 
and the Bechgaard salts, (TMTSF)2X (where TMTSF 
is tetramethyltetraselenafulvalene, and X — PFg, CIO4, 
or Re04), which exhibit simultaneous field-induced spin 
density waves (SDWs) and QHE |,|. 

Electronic transport parallel to the field direction is 
a powerful probe of these chiral surface phases. Our 
analysis reveals that for the integer quantum Hall effect 



(IQHE), the motion along the z-direction is always dif- 
fusive, independent of the impurity scattering strength. 
Localization effects - weak or otherwise - are completely 
absent. The 2d resistivity p^z is temperature indepen- 
dent as r — > 0, with a value which can be much larger 
than h/e^. But perpendicular to the field, along the x- 
axis, the transport is ballistic and the resistivity p^x van- 
ishes as T — > 0. Our predictions for various experimental 
quantities which should be accessible in the Bechgaard 
salts, are given in detail at the end of the paper. 

The surface sheath in the fractional quantum Hall ef- 
fect at filling v — 1/3 is even more anisotropic. In this 
case, we find that electron transport along the z-axis is al- 
ways incoherent. At low temperatures insulating behav- 
ior is predicted with a resistivity diverging as pzz ^ '^/T^, 
even in the absence of any impurity scattering. At T = 
inter-layer transport is completely absent, and the elec- 
trons are "confined" to Id. This chiral surface sheath for 
V = 1/3 is a nice example of a 2d non- Fermi liquid phase. 



FIG. 1. (a): Geometry of a 3d quantum Hall sample. 
2-axis transport is included via the tunneling amplitude t, 
and impurity scattering by the random potential V . (b): As- 
sociated Fermi surface for the sheath. It differs from that of a 
conventional open Fermi-surface metal by the absence of the 
left-moving (dotted) half. 

Consider first the IQHE. The electronic states at the 
edge can be modeled simply as chiral Fermions Q . In the 
3d case of interest here, there is a single edge state per 
layer {x-y plane) per filled Landau level (see Fig. 1). We 
focus on the case v = 1; higher integer filling fractions 
behave similarly. Including a small inter-layer matrix 
element t leads to the non-interacting Hamiltonian 

H = ^ J dx^pjivdxipi - t (^i/jj'ilj.+i + tpj^-^^tpi^ (1) 

i 

= {vpx - 2tC0SPz)lp''{Px,Pz)'fpiPx,Pz), (2) 



where = dpx/{2iT) /^^ dpz/{2n) and V'J is an elec- 
tron creation operator for the edge state in the z*'' layer. 
Here we have set the layer spacing a equal to one. Eq. || 
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shows that inter-layer hopping induces a small disper- 
sion along the layering direction, resulting in "half" of 
an open Fermi surface (see Fig. lb). 

Before including impurity scattering, we re-examine 
the validity of Fermi-liquid theory (FL) with electron- 
electron interactions present. Despite the chiral nature 
of the sheath, the phase-space restrictions which stabi- 
lize the FL for a conventional Fermi-surface continue to 
apply here. In particular, the quasi-particle decay rate 
Iml](a;) ~ w^lnw, the d.c. conductivity a^z ~ l/T'^, 
and the specific heat C j (av) . Long-range Coulomb 
interactions are screened in the static (w — > 0) limit. 

To study d.c. transport at low temperatures, we must 
include impurity scattering at the surface, which we in- 
clude by adding to the Hamiltonian a term of the form. 



Hi, 



E 



dxV{x, iajipjipi 



(3) 



For simplicity, we assume the random potential to be 
Gaussian and uncorrelated, V{x, z)Vj{x' , z') — A6{x — 
2^0^2,2' • contrast to the results for the pure sys- 
tem, the chiral nature of the dirty sheath leads to dra- 
matic differences from ordinary dirty 2d metallic behav- 
ior. In particular, localization along the x-direction is 
precluded a priori since the wavefunctions are nec- 
essarily de-localized along the x-direction. This follows 
from the continuity equation — • J, where 

J(x,z) = {2t/v)lm[<^*{x,z)<^{x,z + a)], valid for arbi- 
trary random potential V . To study possible localization 
effects along the z-axis, we consider the usual averaged 
product of advanced and retarded Green's functions for 
non-interacting electrons, 



V{v;Vt) =G+{Q,r-E)G-{r,0-E + n). 



(4) 



Here G± = [ivd^, - tVl - E t iv + ^(^)]"\ with 
the discrete Laplacian, and rj a positive infinitesimal. A 
standard summation of ladder diagrams, which gives a 
diffusive form for non-chiral Fermions, yields the approx- 
imation 



2-Kp 



-i{n + vpx) + Dpi + 2i]' 



(5) 



where the density of states p = l/27rwa, and D = Pto^, 
with the relaxation time r — 2v/A. The anisotropic form 
of the denominator in Eq. |^ reflects the difference in prop- 
agation along X (ballistic) and z (diffusive). The Einstein 
relation azz — e^pD can be verified by explicit compu- 
tation. The derivation of the diffusive transverse motion 
in Eq. |^ has neglected possibly important quantum in- 
terference corrections. Indeed, for isotropic 2d systems, 
such corrections invalidate diffusive behavior at low tem- 
peratures, due to localization. To study possible interfer- 
ence effects here, we use the Q-matrix approach, which 
allows a systematic treatment of corrections to the sim- 
ple ladder result in Eq. ||. Following McKane and Stone 
1^ , averaged correlation functions are obtained from the 
replicated partition function Z = J [d4i\[d(l)\[dQ] exp(— S*), 
where 



S = Y^J dx^^i(l)A [ivd^ - <Vi -E~ AQ] 



(6) 



puted in the limit n 0. Writing Q 



Here is a 2n-component complex vector, Q is a 2n by 
2n hermitian matrix, and A = 1 cr^. Eq. ^ follows 
from introducing (j)± fields to generate retarded and ad- 
vanced Green's functions, averaging the n^^ moment of 
the original generating function Zq, and decoupling the 
resulting quartic interaction using the Q matrix. For 
r] — 0, this action has a non-compact 0{n, n) symmetry, 
reflecting the physical equivalence of replicas. In many 
respects, rj acts as a symmetry breaking fleld, but is also 
essential to obtain convergence of the functional integral. 
Physical correlators, which are derived from the loga- 
rithm of Zq (InZo ~ lim„^o(^o ~ 1 niust be com- 

?++ Q+ 

Q-+ Q- 

the density of states is, e.g., p = (l/7r)Im(Q++;Qa), 
where no sum is implied over the (arbitrary) replica in- 
dex 1 < a < n. Up to non-singular contact terms, the 
correlator X'(r; 0) = {Q+-.aa{r)Q_+.aa{0))- 

The ladder approximation is recovered as the leading 
term in a saddle point expansion of the effective action 
for Q, obtained by integrating out the bosonic flelds. At 

the mean- fleld level, Q++ — ~Q — i/{2va), reflecting 

the finite density of states. 

Expanding Q — Q + Q, the Q++ and Q fields have 

only massive fiuctuations around the saddle, but the 

and Q ^ fields are massless Goldstone bosons for 77 = 0. 

Following the conventional non-linear sigma model ap- 
proach, we ignore the massive fiuctuations and focus on 
the Goldstone modes, staying within the saddle point 
manifold by eliminating Q++ and Q via the zero mo- 
mentum, tree equations of motion. The resulting theory 
may be expressed solely in terms of the off-diagonal sub- 
matrices of Q, and is governed by an action of the form 

SeS^T^ I i-ivpx + Dpl)Tv[Q+^{p)Q^+i-p)] 
27rp Jp 

+ n / r{p^mQ+-Q-+Q+-Q-+], (7) 

to 0(Qj_^). The interaction vertex ^{Pi) = 
5{Y.Pi)f{Pi)^ with f{p) cipx + C2pI for smaU mo- 
menta, and all higher vertices must have similar momen- 
tum structure owing to the continuous 0{n, n) symme- 
try. Truncating at quadratic order leads to the ladder 
result for P, while retaining the anharmonic corrections 
is analogous to the usual expansion of the appropriate 
non-linear sigma model. 

With this formulation, the stability of the "diffusive" 
metal can be evaluated from simple power counting. 
To make the quadratic action dimensionless, we rescale 
z — > &z, X ^ b'^x, and Q b^^^^Q (with 6 > 1 to 
focus on long- wavelength behavior). Under this trans- 
formation the quartic coupling amplitude Fg ^o/b, 
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implying the irrelevance of anharmonic terms in 

and the stability of the metallic phase. This is in con- 
trast to the usual case of a non-chiral 2d dirty conductor, 
in which the metallic fixed point is marginally unstable. 
Physically, the absence of localization effects may be at- 
tributed to the ballistic motion along the x axis, which 
suppresses multiple scattering. Although it is possible 
that strong enough disorder could induce some kind of 
transverse localization, we have been unable to find evi- 
dence of such a phase, and suspect that metallic behavior 
persists for all impurity concentrations. 

In other 2d electron models which exhibit diffusion 
at T = 0, the de-localized wavefunctions typically ex- 
hibit multi- fractal scaling [||. A classic example is the 
IQHE plateau transition, which manifests multi-fractal 
behavior via anomalous scaling of 'D{p,fl), in the limit 
Dp^ << n. But in the present case 2?(p, $7) scales triv- 
ially with Q and pz , implying that the wavefunctions are 
not multi-fractal. 

Having established the stability of the metallic phase in 
the non-interacting case, we turn to the combined influ- 
ence of interactions and disorder. As shown by Altschulcr 
and Aronov (sec Ref. and Refs. therein), interaction 
effects are much enhanced by diffusive motion in conven- 
tional metals. For example, diffusive relaxation of density 
fluctuations leads to singularities in the tunneling density 
of states. For our semi-ballistic metal, one expects this 
effect to be absent, since any charge buildup is swept 
away at the chiral Fermi velocity. Indeed, one can show 
that the diagrams responsible for the singularity yield 
only smooth contributions in the chiral case. 

The z-axis transport for the surface sheath in the IQHE 
behaves much as a conventional metal. Dramatically dif- 
ferent results are obtained for fractional filling factors, 
as we now describe. We focus on the odd-denominator 
fractional states, particularly ly — 1/3. The edge states 
are then chiral Luttinger liquids of charge v Laughlin 
quasiparticles |Q. Because these quasiparticles have in- 
tegrity only within a single quantum Hall plane, inter- 
layer transport must involve the tunneling of physical 
electrons^ or equivalently l/v quasiparticles. The surface 
"sheath" can be described using the bosonized Euclidean 
action 

/oo ( 1 

dx J dTi^^—dx(f>i{idr(l)t+vdx(l>t) 

-tcos[{(j)i - (8) 

where (/) is a boson field, and /3 ~ 1/T. The charge den- 
sity is ni{x) = dx4'i / {2n) , and the operator e"^'-^'' creates 
a quasiparticle (2tt soliton in (f>) at position x. The co- 
sine term describes inter-layer tunneling with amplitude 
t. In the ideal H \\ z geometry, no flux penetrates be- 
tween successive edge modes; oscillatory phase factors 
are therefore absent from this term. To study the effects 
of tunneling, we employ standard renormalization group 
methods perturbative in t. After introducing a cutoff A 
in Qx, one integrates out modes with Ae~'^^ < \qx\ < A. 



Upon rescaling x —> e x and imaginary time t — > e r, 
the quadratic action is brought back to its original form. 
The tunneling amplitude is, however, renormalizcd 

^==(2-l/.K (9) 

and is technically an irrelevant perturbation. To extract 
the z-axis conductivity imagine renormalizing until the 
rescaled temperature is of order the quantum Hall gap 
Eg, which implies a rescaling factor = Eg/T. From 
scaling one then obtains an inter-layer tunneling time 
T„ ~ T^[Eg/t{t)Y, with t{(.) - tiT/Egf/'^-^. Here 
~ 1/T is a thermal de-phasing time . Since t{t} 
scales to zero at low temperatures. Tin » and the 
inter-layer transport is incoherent. The z-axis conduc- 
tivity can then be obtained from the diffusion constant 
D ~ a?' /Tin and the Einstein relation Uzz = e^pD. This 
gives Gzz ~ {g^ /h){a/vTin) = {const)T'^/'^~^ , which van- 
ishes as T ^ 0. 

We now summarize our predictions for transport and 
thermodynamic quantities. For temperatures well be- 
low the bulk QHE gap, the surface sheath dominates 
the z-axis transport. In this regime, for = 1, a tem- 
perature independent (sheet) conductivity is predicted: 
e^zziv = 1) ^ {e^t^T)/2Tiv. Each additional full Lan- 
dau level will give another surface "sheath" , which con- 
tributes additively to the sheet conductance, so azz{v — 
N) w Nazz{^) (see Fig. 2). This behavior is roughly 
consistent with recent data on (TMTSF)2C104 Hi. As 
discussed above, the v = 1/3 state, if present, has a van- 
ishing z-axis conductivity, Uzz ^ ■ In contrast to the 
quantized Hall plateaus, the "plateaus" in Gzz are un- 
quantized, and will probably exhibit small non- vanishing 
slopes even at the lowest temperatures, as depicted in the 
Figure. However, the "step-like" features, which arise 
from the 3d phase transitions between successive bulk 
Hall phases, can be sharp at T = 0. The low tempera- 
ture behavior of the "steps" depends upon the order of 
the bulk phase transition. In (TMTSF)2X, these transi- 
tions might be first order, with the SDW period chang- 
ing discontinuously. In this case, we expect a discontin- 
uous jump in 17 zz- However, experience with the usual 
QHE leads us to expect continuous critical behavior in 
2d electron gas multi-layers; sufficient disorder may also 
drive the SDW transitions second order. In this case, 
the "steps" will exhibit critical singularities, as we now 
describe. 

FIG. 2. Predicted behavior of azz- We have indicated a 
small with 5B, such as would be present at finite temperature 
in the case of an isotropic bulk critical point. 

A continuous bulk plateau transition is worthy of study 
in its own right. Here we restrict the discussion to 
scaling properties which influence the surface sheath. 
Near the transition one expects an (in-plane) localiza- 
tion length which diverges as ^ ~ \SB\~'^ . Allowing for 
anisotropic scaling, the z-axis localization length is as- 
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sumed to vary as ^ (we expect ( < 1 on physi- 
cal grounds), whereas characteristic energy scales vanish 
as a; ~ with z a dynamical exponent. Following 
standard methods, we find that the bulk conductiv- 
ity crj""* ^'""'^^(T^^)- This form connects activated 
behavior away from the transition (with energy scale 
Tq ^ Eg{iB / £,y , where £b is the magnetic length) to 
the power-law temperature dependence cr^z""^ '-^ 7^(2-0/^ 
at the critical point. As an aside, we point out that 
c^xx ~ T^/^ — > at the 3d critical point, in contrast with 
the prediction of a universal finite value in 2d. 

Upon approaching such a continuous bulk transition, 
the surface sheath develops a thickness ^ ^ ^s, and 
we may estimate the z-axis sheet conductivity as Uzz ^ 
erf Using the bulk scaling results, this implies singu- 
lar behavior at zero temperature for the "steps" in a^z'- 



\SB\ 



(10) 



Experience with other similar disordered systems sug- 
gests that the critical behavior may be isotropic, so that 
C = 1 |11|. If so, we again find a jump in azz- A jump 
due to an isotropic but continuous bulk transition can 
be differentiated from a bulk first order transition by its 
behavior at finite temperature: in the former case, the 
step develops a width SB ^ 7^1/(2")^ while for a first or- 
der transition the discontinuity should be present even 
at r 7^ 0. 

The SDW compounds provide a unique opportunity 
to study the specific heat of quantum Hall edge modes, 
in this case those of the whole sheath. In 2d electron 
gases, the edge contribution to C is masked by that of 
the localized states in the interior. While this should 
also be the case in multi-layers, the SDW gap implies 
that the bulk contribution to C is activated at low tem- 
peratures. One can then measure the specific heat from 
the gapless surface modes. They contribute a linear 
temperature dependence to the specific heat per area: 
Ca — N {irk'^T / Qfiav) . Since the ratio C/T is linear in 
the number N of full Landau levels, one expects "steps" 
in this quantity as well. Recent specific heat data in the 
C/O4 Bechgaard salt appears roughly consistent with this 
variation. The evidence for a latent heat peak near the 
transitions is less clear [T^ ]. 

Finally, we consider the effects of a transverse magnetic 
field {By ^ 0) on the integer surface sheath, induced (al- 
beit non-uniformly) by tilting the sample. In the gauge 
Ax = —ByZ, this field shifts the relative Fermi momenta 
between edge modes in adjacent layers by Skp = eaBy. 
In the absence of impurities, momentum conserving inter- 
layer tunneling is not possible below an energy scale, 
Ab veaBy. This leads to activated z-axis transport, 
Czz ~ Czz.o exp(— As/T), characteristic of decoupled Id 
conductors. With impurity scattering, a finite conductiv- 
ity will be restored at T = 0. But for a relatively clean 
surface, a positive magneto- resistance for z-axis transport 
is expected. 
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